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FOREWORD 


This  report  presents  the  results  of  a  short-term  study 
undertaken  to  fill  a  gap  In  a  rapidly  developing  technology 
of  gliding  parachute  systems.  While  the  problem  is  not  very 
complex,  it  is  hoped  that  a  concise  presentation  of  the  equa¬ 
tions  and  illustrative  calculations  will  he  of  value. 

The  author  wishes  to  express  his  appreciation  to  Mr. 
Richard  J.  Greene  whose  competent  preparation  of  the  illu¬ 
strations  facilitated  timely  completion  of  the  study. 

This  study  w&a  conducted  under  Department  of  the  Army 
Project  No.  1F162203D195 ,  Exploratory  Development  of  Airdrop 
Systems . 
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ABSTRACT 


Equations  are  presented  and  evaluated  for  estimating 
the  wind  effect  on  the  approach  path  and  descending  orbit 
of  gliding  parachute  systems  with  non-proportional  automatic 
homing  control.  Exact  equations  are  presented  for  calculating 
the  overall  downwind  displacement,  the  total  time,  and  the 
positions  of  four  characteristic  points  on  a  typical  des¬ 
cending  orbit,  or  complete  circuit,  over  the  target.  Itera¬ 
tion  equations  are  presented  Incorporating  homing  simulations 
for  calculating  positions  at  equal  time  Intervals  on  the  In¬ 
itial  approach  and  descending  orbit.  The  control  response 
time,  effect  of  variations  In  deployment  position,  and  Impact 
probability  are  discussed.  The  equations  are  evaluated  for 
systems  having  horizontal  velocities  of  40  and  60  fps  and 
having  turn  radii  of  75  and  100  ft  in  wind  speeds  of  from 
4  to  32  fps . 

Based  on  results  of  a  computer  study,  empirical  equa¬ 
tions  are  presented  for  calculating  the  probability  of 
Impact  within  a  given  radius  of  the  target  and  for  cal¬ 
culating  the  radius  of  the  circle  of  equal  probability  as 
functions  of  turning  radius,  wind  velocity,  and  horizontal 
system  velocity  relative  to  air  (or  glide  ratio  assuming 
standard  airdrop  descent  velocity)  .  The  analysis  shows  that 
Increasing  horizontal  system  velocity  relative  to  air  decreases 
wind  distortion  both  of  a  crosswind  approach  and  of  the 
descending  orbit  thereby  Increasing  the  probability  of  Im¬ 
pact  close  to  the  target. 


viii 


la  Introduction 
a.  Purpose 


Tha  practical  design  of  gliding  parachute  systems  ra- 
qulraa  estimation  of  performance  in  a  wind  whlla  tha  system 
la  unuair  autoaatic-hoaint  radio  control.  Vhlla  aaat  gliding 
ayatoaa  would  hava  a  aanual  control  capability  which  would 
allow  a  nan  on  the  ground  to  correct  for  wind  of facta ,  auto¬ 
matic  control  ahould  ba  aufflclantly  accurata  to  facilitate 
una  by  relatively  untralnad  and  inexperienced  flald  peraonnel 
and  to  facilitate  uaa  in  Halted  visibility.  Moat  gliding 
parachute  ayataaa  hava  glida  ratios  froa  2.0  to  3.0  end  hava 
a  rata  of  daacent  of  about  20  ft/aae-  Tharafora,  tha  hori¬ 
zontal  velocity  eoaponant  ia  generally  betwaen  40  and  60 
ft/aae.  A  typical  wind  speed  of  13  knots  or  25  ft/sae  Would 
have  a  significant  effect  on  the  ground  track  of  such  sys¬ 
tems.  This  report  presents  both  basic  equations  for  deter¬ 
mining  certain  ground  track  characteristics  and  an  extensive 
computar  study  which  shows  tha  quantitative  effects  of  various 
systam  parameters  on  tha  ground  track  in  a  wind.  Also,  an 
empirical  aquation  la  presented  for  uaa  in  design  of  tha  con¬ 
trol  systam  to  satisfy  given  performance  requirements . 

b.  Scope  of  the  Analysis 

The  analysis  is  based  on  a  non -proportional  radio  con¬ 
trol  system  which  would  Initiate  control  only  for  a  full 
left  turn,  full  right  turi,  or  straight  flight  dapending 
on  orientation  of  the  system  toward  tha  transmitter  which 
is  located  at  the  target  point.  This  systam  is  lass  expen¬ 
sive  and  mechanically  simpler  than  a  proportional  system. 

Also,  the  author  is  awera  of  ona  control  systam  of  this  type 
currently  in  experimental  usa. 

Equations  are  prasantad  in  this  paper  both  for  cal¬ 
culating  certain  characteristic  points  on  tha  final  des¬ 
cending  orbit  and  for  computing  Instantaneous  points  on  the 
approach.  Simple  analyses  of  velocity  vactors  end  geometry 
ara  used  to  calculate  points  on  a  descending  orbit  and  the 
time  for  ona  comp  lata  orbit,  or  circuit,  ovar  the  target. 

By  considering  small,  time  increments  during  controlled 
turns,  tha  ground  track  during  homing  approach  and  descending 
orbits  can  ba  approximate.  A  basic  estlmata  of  tha  deploy¬ 
ment  ahvelope  for  an  offsat  delivery  mission  is  presentad  to 
show  tha  range  of  possible  error  in  altitude  of  the  system 
at  the  completion  of  the  initial  approach  to  the  targat. 

Basad  on  the  assumption  that  Impact  will  occur  any  time  during 
a  typical  descending  orbit  ovar  tha  target,  impact  position 
probability  is  considered  aquel  to  the  time  spent  within  a 
glvan  distance  of  the  targat  in  percent  total  time  for  a 
complate  orbit. 


2.  Methods  »f  Analysis 
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Effects  of  wind  on  the  performance  of  gliding  systems 
must  be  snelyzed  from  the  concept  that  the  system  has  a  con¬ 
stant  total  velocity  relative  to  tha  air  regardless  of  wind. 

When  the  system  executes  e  constant-rate  turn  for  a  given 
amount  of  time,  the  change  in  its  position  over  the  ground 
may  be  determined  by  calculating  the  position  it  would  have 
without  wind  end  then  displacing  the  position  downwind  a 
distance  equal  to  the  product  of  the  steady  wind  speed  and 
the  elapsed  time.  Orientation  of  the  system  1s  parallel 
to  tho  orientation  it  would  have  without  wind. 

a.  Beale  Orbit  Geometry 

Four  characteristic  points  and  the  total  time  of  the 
typical  descending  orbit  can  be  calculated  from  exact  re¬ 
lations.  The  results  apply  accurately  to  the  second  and 
subsequent  passages  over  the  target  which  are  Independent 
of  the  initial  approach.  Consider  the  geometric  relations 
shown  in  Figure  1  where  V  is  the  horizontal  velocity  com¬ 
ponent  of  the  system  relative  to  air*  W  is  the  wind  speed 
(parallel  to  the  Y-axis)  and  r  is  the  radius  of  the  turn 
of  the  system  in  calm  air.  Circular  motion  to  the  left  of  < 

the  Y-axis  is  assumed.  The  positions  relative  to  the  X- 
axls  end  the  time  relationships  of  the  four  points  on  the 
actual  orbit  are  derived  from  the  positions  on  an  analogous 
circular  orbit. 


Fig  1.  BASIC  GEOMETRIC  RELATIONS  FOR  DESCENDING  ORBIT 

? 


Point  A  coincides  with  the  target.  On  a  typical  orbit, 
the  system  would  head  upwind  at  Point  A  and  hold  a  constant 
turn  control  through  Points  B,C,D  and  E  to  return  to  Point  A. 
At  PointB  :,  the  vector  sum  of  the  system  velocity  and  the 
wind  velocity  yields  a  resultant  having  pure  crosswind  direc¬ 
tion.  Hence 

sinOC  *  W/V 

or  OL  ■  arcsin  W/V 

where  Ot  is  the  angle  of  inclination  of  the  vector  V  from 
the  x-axis.  From  the  geometry  related  to  a  ci-cle,  we  see 
that  Oc  is  also  the  angle  from  the  diameter  parallel  to  the 
Y-axis  of  a  line  passing  through  the  point  on  the  circle  anal¬ 
ogous  to  Point  B.  The  angularity  of  the  vector  V  thus  changes 
by  the  amount  7T/ 2  -(^between  Point  A  and  Point  B.  If  Tq  is 
the  time  required  to  turn  a  complete  circle  in  calm  air,  then 
T0  (iT/2  -oc  )  /2ft  is  the  time  required  to  move  from  Point  A 
to  Point  B.  The  coordinates  of  Point  B  may  thus  be  calculated 


from  the  relations 

XB  -  -r  (1  -  sin  Ot  )  (3) 

Yb  -  r  cost*  -  WTC  (  TT/2  -  Ot  )  / 2TT  (4) 

At  Point  C,  motion  is  purely  downwind.  The  coordinates 

of  Point  C  are 

Xc  -  -2r  (5) 

Yc  -  -WT0/2  (6) 

where  the  system  has  cojplated  a  180-degree  turn.  At  Point 
D,  the  system  has  an  upwind  velocity  component  sufficient  to 
cancel  the  wind  velocity.  The  coordinates  of  Point  D  are 

XD  -  XB  (7) 

Yd  -  -r  cosOC  -  WT0  (3TT/2  +  OC  )  /2TT  (8) 

At  Point  E,  motion,  is  again,  purely  upwind  toward  the  target. 
The  x  coordinate  of  Point  E  is  zero,  and  the  y  coordinate  is 

yE  -  -wt0  (9) 

The  total  distance  from  Point  B  to  Point  D  is 

L  -  2r  cosOt  +  WT0  (  IT  +  2CC  )  /2 IT  (10) 

The  total  elapsed  time  for  the  complete  orbit  is 

T  -  Tc  +  fyEi  /  (V-W)  (11) 


(1) 

(2) 
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or,  substituting  the  value  of  yj?  from  Eq  (9), 

T  -  T0  +  WT0  /  n- W)  (12) 

The  value  of  T0  may  be  calculated  from 

T0  -  2  'TTr /V  (13) 

This  value  may  be  substituted  into  each  of  the  above  equa¬ 
tions  involving  TQ  to  give  the  following  simplified  relations: 

yB  -  r  co soc  -  r  (W/V)  ( -TC  /2  -0C  )  (14) 

yc  --TTr  (W/V)  (15) 

yD  -  -r  cosOC  -  r  (W/V)  (37T/2  +0()  (16) 

yE  -  -2  TT  r  (W/V)  (17) 

L  -  2r  cosOC  +  r  (W/V)  ( +  20C  )  (18) 

and 

T  -  2 TTr / (V-W)  (19) 


Note  that,  in  each  geometric  relation,  the  sole  velocity 
term  is  the  dimensionless  ratio  W/V  It  is  interesting 
that  the  time  for  a  complete  orbit  in  a  wind  is  equal  to 
the  time  required  to  turn  a  perfect  circle  at  the  minimum 
ground  speed,  V-W.  For  gliding  parachute  systems,  values 
of  W/V  may  commonly  be  of  the  magnitude  of  1/2.  Hence,  the 
above  equations  show  that  significant  distortions  of  the 
calm  air  orbit  may  commonly  occur. 

The  four  characteristic  points,  the  overall  length,  and 
the  orbit  time  are  sufficient  for  basic  estimation  of  wind 
effects.  However,  for  any  time,  t,  during  an  orbit,  the 
position  of  the  system  can  be  determined  from  the  following 
equations : 


For 

Oit£T0, 

x  * 

-r  (1  -  cos  (Vt/r)) 

(20) 

y  - 

r  sin  (Vt/r)  -  Wt 

(21) 

For 

T0  <  t  T ,  x  ■  0  and 

y  - 

-WTq  +  (t  -  T0)  (V-W) 

(22) 

However,  for  calculating  instantaneous  positions,  the  equa¬ 
tions  presented  below  in  Section  2b  are  more  suitable  since 
they  apply  to  the  approach  as  well  as  to  the  descending  orbit. 
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From  the  basic  relations,  we  see  that  distortion  in¬ 
creases  In  direct  proportion  to  r  and  In  inverse  proportion 
to  V  and,  hence,  to  the  glide  ratio.  Maximum  glide  ratio 
would,  therefore,  minimize  distortion  for  a  given  radius  of 
turn . 


Fig  2.  GEOMETRY  OF  TURN  DURING  APPROACH 
(WITHOUT  WIND  DISPLACEMENT) 


b.  Iteration  Equationa  for  Homing 

Points  at  small  time  intervals  along  the  path  of  the 
homing  approach  may  be  determined  from  the  assumption  that, 
at  any  position,  the  system  can  turn  through  a  small  are  of 
angle  0  either  to  the  left  or  to  the  right  -  whichever  brings 
It  closer  to  the  target.  Geometric  relations  are  derived 
first  for  the  case  of  calm  air.  The  wind  effect  Is  added 
as  a  constant  displacement  at  the  end  of  each  arc.  The 
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tangential  velocity  Is  equal  to  the  constant  horizontal 
velocity  of  the  system  since  the  system  has  aerodynamic 
equilibrium.  Actually,  the  velocity  and  glide  ratio  may 
be  less  in  a  turn  than  in  straight  flight.  However,  it 
will  be  shown  later  that  straight  flight  is  virtually  im¬ 
possible  for  the  assumed  control  system.  Hence,  the  velo¬ 
city  may  be  considered  constant.  From  the  geometry  shown 
in  Fig  2,  the  position  and  orientation  of  the  system  upon 
completion  of  the  arc  9  may  be  calculated  exactly. 

Assume  the  radius  of  turn,  r,  the  arc  angle,  9,  and 
the  horizontal  velocity  component,  V,  relative  to  air  remain 
constant.  Given  an  initial  position,  (x,y),  and  initial 
velocity  components,  Vx  and  Vy ,  which  indicate  orientation 
relative  to  the  target  at  the  origin  of  the  coordinate  system, 
the  two  possible  end  points  (xi,  y^)  and  (x2  ,  y2)  after  com¬ 
pleting  a  left  or  right  turn  through  angle,  9,  may  be  cal¬ 


culated  from  the  relations 

x  -  x^  ■  2r  sin  (9/2)  cos  (0  -  9/2)  (23) 

-(y  -  yi>  »  2 r  sin  (9/2)  sin  (0  -  9/2)  (24) 

x  -  x2  -  2r  sin  (9/2)  cos  (0  +  9/2)  (25) 

-(y  -  y2)  *  2r  sin  (9/2)  sin  (0  +  9/2)  (26) 

where  0  is  the  initial  angle  of  inclination  of  V  from  the 
x-axis.  The  velocity  components  at  points  (xj.,  y^)  and 
(x2 *  y2)  may  be  calculated  from 

Vy  -  V  sin  (0  -  9)  (27) 

V„  -  V  cos  (0  -  9)  (28) 

1 

Vv  -  V  sin  (0  +  9)  (29) 

y  2 

VX2  -  V  cos  (0  +  9)  (30) 

The  angle  of  0  has  the  relations 

sin  0  -  Vy/V  (31) 

cos  0  =  -Vx/V  (32) 

Where  Vy  and  Vx  are  the  values  prior  to  the  turn. 


If  Si  and  S2  are  the  squares  of  the  distances  from  the 
origin  to  points  (x^  ,  yi)  and  (X2  ,  y2)  as  defined  by 
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«n»: 


S1 

s„ 


‘1 


?l‘ 


x22  +  y  2  2 


(33) 

(34) 


then  homing  control  is  effected  by  choosing  the  shorter  dis¬ 
tance  on  the  basis  of  the  sign  of  S2  -  S,.  After  alge¬ 
braic  manipulation  of  Eq  (23)  ,  (24)  ,  (25)  and  (26)  we  find 


Thus  if 

xVy  -  yVx  >0 

than 

S2  -  Sx  >  0  ; 

and  if 

xVy  -  yVx  < 

than 

S  -  S.  —  0  . 
2  1 

In  final  form,  the  equations  for  the  coordinates  x' 
and  y'  of  the  end  point  of  the  arc  and  the  corresponding 
velocity  components  Vx'  and  V  '  ire  given  by 


X 

xS  -  x  + 

(2r/V)  sin  (9/2)  [Vx  cos  -  aVy 

sin  |] 

(35) 

y  '  m  y  + 

—•  sin  |  [Vy  cos  ^  +  aVx  sin  p 

-  wr 

(36) 

• 

Vx'  -  V 
x  X 

-  2  sin  |  (Vx  sin  |  +  aVy  cos  p 

(37) 

> 

1 

> 

+  2  sin  5.  (aV  cos  %  -  Vv  sin  %) 
2  x  2  y  * 

(38) 

+1  if  (xVy  -  yVx)  ^  0 


where  a 

y  •  • 

and  a  -  -1  if  (xVy  -  yV  )  —  0. 


These  equations  may  be  used  in  an  iterative  calculation 
to  find  all  subsequent  points  on  the  ground  track  during 
the  approach  and  the  final  descending  orbit.  The  sole 
correction  for  wind  effects  is  made  by  subtracting  the 
term  WT*  from  the  y  coordinate  as  shown  in  Eq  (36)  where  T 
is  the  time  elapsed  during  the  turn.  Values  of  Vx'  and 
Vy'  are  relative  to  air  and  remain  unaffected  by  wind. 
Hence,  the  values  of  x',  Vx',  Vy  '  and  the  corrected  value 
of  y'  may  be  used  to  determine  the  value  of  a  for  the 
subsequent  turn. 

The  angle  0/2  may  be  calculated  in  terms  of  the  time 
increment  by  the  relation 


<?> 


180  V  y/2  TTr 


(39) 


-  -1 


for  9/2  In  degrees. 

c.  Significance  of  the  Iterating  Increment 

In  specific  applications,  the  value  of *X  used  in  Eq  (39) 
may  be  chosen  to  represent  a  characteristic  response  time  of 
the  parachute  system.  Prior  to  a  turn,  the  physical  parachute 
system  must  undergo  a  lateral  oscillation  governed  by  system 
geometry  and  by  the  effective  masses  of  the  canopy  and  of 
the  suspended  load.  While  a  detailed  analysis  of  this  phen¬ 
omenon  is  beyond  the  scope  of  this  report,  general  consider¬ 
ations  show  that  the  period  of  the  oscillation  is  a  function 
of  parameters  which  would  remain  virtually  constant  through¬ 
out  gliding  flight.  The  analysis  presented  in  the  report 
"Dynamic  Stability  of  a  Parachute  Point-mass  Load  System" 
by  H.G.  Heinrich  and  L.W.  Rust  Jr.(^  shows  comparable  dy¬ 
namic  relations.  A  similar  analysis  could  be  used  to  estimate 
the  response  time  of  a  gliding  system  to  a  turn  control.  If 
the  response  time  for  a  particular  system  can  be  determined, 
the  use  of  such  a  value  for  the  iterating  increment  of  time, 

»  *n  ^9  (39)  would  increase  the  accuracy  of  the  computed 
ground  track.  The  physical  system  would  be  capable  of  ex¬ 
ecuting  homing  control  only  at  equally  spaced  time  intervals 
according  to  the  oscillatory  motion. 

For  the  computer  evaluation  presented  in  this  report, 
the  values  of'T  range  from  0.3  to  2.4  sec.  This  range  seems 
reasonable  in  view  of  the  known  characteristics  of  standard 
vertically-descending  parachute  systems.  For  example  the 
time  required  for  a  system  consisting  of  a  5,000  lb.  payload 
and  100-ft.  diameter  circular  canony  to  rotate  from  27°  to 
0°  vertical  inclination  is  approximately  2.3  sec  as  shown 
by  flight  test  data.  Personnel  parachutes  exhibit  similar 
oscillation  times.  However,  the  values  of  T  used  in  this 
computer  study  were  calculated  as  4 ,  8  and  16%  T  to  give 
incremental  turns  of  14.5,  29  and  58  degrees.  Values  of  Y 
less  than  4%  T0  required  impractical  computation  times  and 
did  not  result  in  significantly  smoother  homing  patterns  as 
discussed  in  section  3a.  For  values  of'T  greater  than  16% 

T0  and  approaching  25%  Tc,  the  computational  methods  of  sec¬ 
tion  2b  show  adverse  peculiarities.  The  assumption  that  the 
track  remains  circular  during  the  time'Tis  probably  Invalid 
for  larger  values  of'T  ,  since  the  actual  path  is  somewhere 
between  a  straight  tangential  line  and  a  pure  circular  arc 
during  the  oscillatory  motion. 


^^Report  No  FDL-TDR-6  4- 1  n  ''June  1965)  USAF  Flight  Dynamics 
Laboratory,  WPAFB ,  Ohio. 
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d.  Effect  of  Deployment  Envelope 


The  range  of  variation  in  altitude  of  the  gliding  system 
as  it  first  arrives  over  the  target  shows  that,  in  general, 
several  orbits  or  360-degree  turns  will  be  completed  before 
impact.  Obviously,  if  the  deployment  position  could  be  con¬ 
trolled  sufficiently,  ground  impact  would  occur  at  the  target 
regardless  of  the  orbit  dimensions  since  final  orbiting  could 
be  avoided.  Although  the  glide  ratio,  X  /d ,  in  aerodynamic 
equilibrium  is  constant,  random  wind  variations  and  random 
system  orientation  during  deployment  must  be  assumed  for  a 
mission  which  utilizes  the  full  offset  capability  of  the  gliding 
system.  Surface  winds  will  seldom  be  constant  between  the 
launch  and  target  positions.  Let  A  W  be  the  maximum  random 
wind  velocity  which  may  occur  either  toward,  or  away  from, 
the  target.  Then  the  effective  glide  ratio  has  the  range  of 
from 


(  J^/d) 


V  -AW 
Vd 


(^/d) 


V  +AW 
Vd 


where  Vd  is  the  constant  rate  of  descent  and  V  is  the  constant 
horizontal  velocity  relative  to  air.  Navigational  error, 
varying  loss  of  altitude  during  deployment,  and  random  orien¬ 
tation  of  the  system  at  completion  of  deployment  will  contribute 
to  uncertainty  in  range  and  altitude  at  the  instant  the  system 
starts  gliding  directly  toward  the  target.  While  a  complete 
analysis  of  these  factors  is  beyond  the  scope  of  this  report, 
assume  the  range  is  i  0  +  and  hQ  +  Ah  is  the  altitude  when 
the  system  is  fully  deployed  and  initially  oriented  toward  the 
target.  The  variation  in  altitude  at  the  first  pass  over  the 
target  will  then  be  approximately  between 

hmin  "  (ho  “  Ah)  -  '  S.  .  *  AJb  ' _ Vd _ )  (A2) 

V  -  A  W 


Ah)  - 


Thi3  assumes  that  launch  is  at  a  prevailing  upwind  location 
since  a  crosswind  approach  would  be  subject  to  additional 
variations.  To  plan  a  typical  offset  mission  for  the  gliding 
system,  a  value  of  Q  would  probably  be  chosen  on  the  basis 
of  tactical  requirements.  Values  of  Ah,  &X  ,  V<j ,  V  and  aW 


would  be  known  or  estimated.  A  i 
found  from  Eq  (42)  assuming  hrajn 
value  of  h _  is 


suitable  value  of 


could  be 


Then  the  corresponding 
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(44) 


^max 


2  Ah  +  2  Vd 


'  +  vA*  ) 

(V  +  A  W)  (V  -  A  W) 


and  the  median  altitude  is 

h  -  ah  ♦  vd  (45) 

(V  +  A  W)  (V  -  A  W) 

To  determine  the  order  of  magnitude  of  h  consider  the  following 
s ample  values : 

0  -  35,000  ft. 

=  500  ft. 

%h°  =  300  ft. 

V  =  60  fps 
Vj  =  20  fps 
AU  *  20  fps  and  0 

For  AW  ■  20  fps,  Eq  (45)  gives  the  value  h  *  4,862  ft  with  hQ  = 
17,550  ft  from  Eq  (42).  If  we  assume  the  radius  of  turn,  r,  is 
100  ft,  Eq  (19)  gives  the  value  T  *  15.6  sec.  With  V<j  ■  20.  ft. 
the  altitude  loss  per  orbit  is  312  ft.  The  high  value  of  h  is 
mostly  a  function  of  AW  since  h  «  467  ft  for  -  0.  Thus,  we 

see  that  an  average  of  from  1.5  to  15.6  complete  orbits  would  be 
executed  prior  to  impact  for  A*.  W  between  0  and  20  fps.  There¬ 
fore,  impact  accuracy  is  largely  dependent  on  the  orbit  geometry. 


e.  Impact  Position  Probability 

The  analysis  of  the  previous  section  shows  that,  in  gener¬ 
al,  impact  may  occur  at  any  time  during  a  typical  orbit  over 
the  target.  In  other  words,  impact  will  occur  at  equal  proba¬ 
bility  with  time  during  an  orbit.  Hence,  the  probability  that 
the  system  will  impact  within  a  given  radius  of  the  target  is 
equal  to  a  ratio  of  time  spent  within  the  radius  to  the  total 
time  elapsed  for  a  typical  orbit.  For  the  computer  studv 
presented  in  Section  3,  impact  probability  is  calculated  from 
the  number  of  points,  indicating  equal  time  increments,  hich 
are  found  within  a  given  radius  of  the  transmitter  (i.e.  - 
the  origin  of  coordinates)  as  a  percent  of  the  total  number 
of  points  for  a  complete  orbit.  By  plotting  the  probability 
as  a  function  of  radius  from  the  target,  an  estimate  of  the 
radius  of  the  circle  of  equal  probability  (CEP)  is  obtained. 

An  attemot  was  made,  using  the  equations  of  Section  2a,  to 
derive  an  approximate  equation  for  the  radius  of  the  CEP  since 
the  parameter  is  commonly  used  to  indicate  impact  accuracy 
and  is  of  prime  importance  in  determining  the  size  of  the  drop 
zone  required  by  a  gliding  system.  No  practical  approximation 
was  found  from  these  equations.  However,  the  recurring  appear¬ 
ance  of  the  dimensionless  terms  Rc/r  and  W/V  -  where  Rc  denotes 
radius  of  the  CEP  -  prompted  an  empirical  analysis  of  the  com¬ 
puter  results  in  relation  to  these  parameters.  A  suitable  em¬ 
pirical  equation  was  thus  obt-Hn^'.  Techniques  and  results  are 
presented  fully  in  Section  3d. 

1° 


a 
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3. 


Results  of 


Computer 


Study 


A  computer  study  was  made  of  the  relations  presented  in 
Section  2b  (Fig  3).  The  ground  track  was  computed  from  the 
initial  point  until  completion  of  two  orbits  or  until  a  complete 
orbit  was  clearly  defined.  Although  actual  impact  position 
would  depend  on  the  altitude  at  the  initial  point,  the  system 
would  follow  the  computed  ground  track  prior  to  impact  regard¬ 
less  of  altitude  providing  glide  ratio  and  wind  speed  did  not 
change  due  to  altitude  effects.  In  all  cases  the  wind  was 
assumed  constant  at  the  specified  speed  directed  toward  the 
target  along  the  positive  y  -  axis. 

Comparison  was  made  between  the  iteration  results  and  the 
geometric  relations  presented  in  Section  2a.  Impact  proba¬ 
bility  was  estimated  according  to  the  hypothesis  discussed 
in  Section  2e . 


a .  Cases  Studied 


The  cases  studied,  as  described  in  Table  1,  may  be  con¬ 
sidered  typical  of  most  gliding  parachute  systems.  A  wind 
speed  of  15  knots,  or  25  fps,  is  assumed  for  Cases  1  -  18  as 
a  standard  maximum  operating  condition.  Since  cargo  safety 
during  impact  generally  requires  a  rate  of  descent  of  20  fps, 
horizontal  velocities  relative  to  air  of  40  to  60  fps  were 
chosen  to  exemplify  the  common  range  of  glide  ratio  from  2.0 
to  3.0.  Values  of  the  radius  of  turn  were  chosen  as  75  and 
100  ft.  which  are  reasonable  in  terms  of  aerodynamic  perfor¬ 
mance  and  practical  accuracy  requirements.  However,  gliding 
systems  for  payloads  beyond  1000  lb  will  probably  have  larger 
turn  radii  since  the  turn  radius  is  somewhat  related  to  the 
distance  between  the  canopy  and  the  payload.  Values  for  the 
time  lag ,  T' ,  or  iterating  time1 increment ,  were  chosen  as  4%, 
8%  and  16%  of  the  total  time,  T0 ,  to  turn  360  degrees  in  calm 
air.  In  one  case,  not  included  in  the  final  presentation, 
a  value  of  "T  •  0.4%  of  T0  was  used.  However,  the  results 
differed  negligibly  from  those  for'T-  4%  T0 ,  and  computation 
time  was  extremely  long.  For  cases  1-9,  comparison  is  made 
between  approaches  near  the  target  from  upwind,  crosswind, 
and  downwind  positions.  All  remaining  cases  show  either  the 
approach  from  upwind  or  the  orbit  from  an  initial  position 
over  the  t  arget  . 
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b .  Observations 


The  computed  ground  track  for  the  calm  air  condition  is 
shown  in  Figure  4.  Figures  5  to  23  snow  the  effects  of  wind 
for  Cases  1  to  IS  as  listed  in  Table  1.  Figures  5,8  and  12 
show  that  an  approach  from  an  initial  position  upwind  of  the 
target  causes  the  system  to  execute  a  small  orbit  between  the 
first  and  second  passes  over  the  target.  In  some  figures  an 
oscillating  curve  has  been  drawn  through  points  which  seem 
to  lie  on  a  stright  line.  The  slope  of  the  oscillating  curve 
coincides  with  the  system  orientation  at  each  point  as  deter¬ 
mined  from  the  velocity  vectors  Vx  and  Vy  .  This  is  indicative 
of  the  type  of  motion  although  the  exact  amplitude  is  not  known, 
indeed,  oscillating  orientation  is  not  necessarily  indicative 
of  an  oscillatory  path  over  the  ground.  The  figures  where 
nr  •  SI  and  16%  T0  show  some  pecularities  resulting  from  the 
discrepancy  between  actual  system  orientation  and  the  apparent 
direction  of  the  ground  track.  Figure  11  shows  the  actual 
system  orientations  indicated  by  the  resultant  of  Vx  and  Vy 
for  an  expanded  portion  of  the  ground  track  of  Case  6.  The 
portion  shown  in  Figure  11  occurs  during  the  second  approach 
to  the  target  shown  in  Figure  10.  The  longer  time  lags  in 
several  instances  cause  s  lort  oscillations  at  close  proximity 
to  the  target  as  exemplified  in  Figures  11,12,14  and  20.  Indeed, 
a  peculiar  resonance  is  shown  in  Figures  14  and  20  which  would 
greatly  Increase  the  probable  accuracy.  This  phenomenon  would 
seem,  however,  difficult  to  achieve  and  control  in  practice. 

In  each  of  the  figures  ,  a  circle  of  200  ft  radius  Is  shown  for 
comparison  to  the  calm  air  accuracy  with  r  ^  100  ft. 

Figures  24,25  and  26  show  the  effect  of  wind  on  long 
crosswind  approaches.  For  various  offset  distances  up  to 
2,000  ft  with  the  system  values  of  Case  8,  the  initial  aoproach 
to  the  target  is  consistently  from  the  five-o'clock  direction. 
Maximum  downwind  displacement  increases  from  the  value  for  a 
500  ft  offset  at  a  rate  of  15  ft  per  100  ft  additional  offset 
(Fig  24).  This  rate  is  proportional  to  W/V,  however.  The 
path  for  a  small  time  lag  (Fig  25)  is  virtually  the  mean  of 
the  path  for  the  larger  time  lag. 

Downwind  drift  on  the  approach  is  much  more  sensitive 
to  glide  ratio  than  to  radius  of  turn  assuming  a  constant  rate 
of  descent.  Since  most  gliding  systems  experience  a  loss  of 
glide  ratio  as  radius  of  turn  is  decreased,  the  important 
feature  of  Figure  26  is  that  a  horizontal  velocity  with  large 
radius  of  turn  may  yield  less  distortion  than  a  lower  hor¬ 
izontal  velocity  with  smaller  radius  of  turn.  The  large  time 
lags  were  used  in  Figures  24  and  26  purely  for  shorter  compu¬ 
tation  times. 
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c. 


Correlation  with  Basic  Orbit  Geometry 


The  values  of  W,V,  and  r  for  Cases  1  to  18  were  used 
In  the  equations  of  Section  2a  to  calculate  the  maximum  down¬ 
wind  displacement  (L) ,  the  total  elapsed  time  fT)  oer  orbit 
and  the  coordinates  of  the  four  characteristic  points  (B,C, 

D,  and  E)  for  typical  orbits  corresponding  to  the  computed 
ground  tracks.  The  calculated  results  are  compared  to  the 
measured  values  from  Cases  1,2,3,10,13  and  16  in  Tables  2A 
and  2B.  The  actual  values  as  determined  from  Figures  6,7,8, 
16,10  and  22  are  shown  in  Table  2A  for  comparison  to  the  cal¬ 
culated  values.  The  measured  values  of  T.  and  T  have  good 
correlation  with  the  calculated  values  while  the  other  values 
appear  to  have  poor  correlation  in  several  cases.  However, 

Table  2B  shows  that  the  apparent  position  discrepancy  Is  the 
result  of  delays  in  the  start  of  the  orbits  due  to  inefficient 
homing  commands  when  the  system  passes  over  the  target.  The 
values  for  each  case  shown  in  Table  2B  were  obtained  from  the 
corresponding  values  in  Table  2A  by  adding  constant  displace¬ 
ments  to  the  point  coordinates  so  that  Point  B  was  made  coin¬ 
cident  with  the  calculated  position.  The  resulting  correlation 
of  Points  C,  D  and  E  is  very  good. 

When  the  svstem  returns  to  the  target,  it  generally  passes 
over  the  target  during  one  time  Increment.  Hence,  the  orien¬ 
tation  of  the  system  after  passing  the  target  is  dependent 
on  the  homing  control  exerted  at  the  beginning  of  the  time 
interval  before  passing  the  target.  Although  the  theoretical 
orbit  is  the  result  of  a  constant  rate  turn  in  one  direction, 
the  svstem  orientation  at  -he  end  of  a  time  interval  may 
dictate  a  homing  turn  in  the  opposite  direction7 because  of 
close  proximity  to  the  target.  The  inconsistency,  therefore, 
of  one  or  two  turns  in  opposite  directions  results  in  an  up¬ 
wind  displacement  of  the  point  at  which  the  system  commences 
a  constant  rate  360-degree  turn. 

d,  Effect  of  Parameters  on  Impact  Position  Probability 

The  effect  of  approach  direction  and  time  lag  on  Impact 
accuracy  related  to  a  200  ft  radius  circle  for  Cases  1-9  is 
shown  in  Table  3.  The  important  observations  from  Table  3  are 
that  probability  during  the  first  orbit  is  best  from  an  upwind 
approach  and  that  probability  generally  increases  with  time 
lag.  Figures  12  and  14  illustrate  the  behavior  near  the  target 
accounting  for  the  listed  probabilities  of  100%. 

For  Cases  1,10,13  and  16  in  which  HT  -  4%  T0 ,  probabilities 
were  determined  for  target,  circles  of  50,100,150  ,200,250,300  and 
350  ft  radius.  The  radii  of  the  circles  of  equal  probability 
(i.e.  -  50%  probability)  were  determined  from  a  graph  of  pro¬ 
bability  vs  radius  from  target.  The  results  showed  peculiar 
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relations  and  a  high  degree  of  linearity  which  warranted  fur¬ 
ther  study.  Hence,  three  more  cases.  Cases  19,  20  and  21  of 
Table  1,  were  computed  to  extend  the  range  of  wind  velocity 
studied.  Figure  28  shows  the  resultin'1  ground  tracks.  The 
results  of  the  subsequent  analysis  are  shown  in  Table  4.  Since 
high  dependence  on  the  dimensionless  values  R/ r  and  W/V  was 
anticipated,  the  results  shown  in  Table  4  were  retabulated  in 
Table  4  according  to  R/r  and  W/V.  The  results  from  Table  5 
are  plotted  in  Figure  29  with  corresponding  least  squares  lin¬ 
earisations  at  each  value  of  W/V.  To  obtain  an  empirical  equa¬ 
tion  relating  probability  (P)  to  R/r  and  W/V,  a  linear  least 
squares  analysis  was  made  of  the  variation  with  W/V  of  the 
slopes  (m)  and  of  the  intercepts  (b)  for  each  P  vs  R/r  linear¬ 
ization.  The  analysis  is  illustrated  in  Figure  30  which  shows 
only  limited  conformity  to  the  linear  results.  However,  as 
a  first  approximation,  the  linear  snunre  -  method  may  be 

considered  sufficient.  The  resulting  empirical  function  is 

P  -  (-35.477  £  +  45.085)  |  +  (29.735  H.  -  10.374)  (46) 

or 

P  -  45.085  |  -  35.477  *  ^  +  29.735  |  -  10.374  (47) 

where  P  is  the  percent  probability,  R  is  the  radius  from  the 
target,  r  is  the  turning  radius,  W  is  the  wind  speed  and  V  is 
the  horizontal  system  velocity  relative  to  air.  If  we  let 
P  -  50%,  we  can  solve  Eq  (47)  for  the  radius  of  the  CEP,  Rc : 


Rp/r 


60.374  -  29.735  W/V 
45.085  -  35.477  W/V 


Rc/r 


1.70  -  0.84  W/V 

1727  -  w/v - 


The  curves  from  Eq  (47)  superimposed  on  the  data  in  r‘ guro  r,l 
show  particularly  good  agreement  at  p  =  50%.  Therefore,  the 
linear  approximations  yield  sufficient  accuracy  as  resolved  in 
Eq  (49).  A  plot  of  Eq  (49)  is  presented  in  Figure  1 1  )  0  W/V£ 

1.0.  Good  agreement  with  theory  is  found  at  the  end  point, 

W/V  »  0  and  W/V  »  1.0.  For  W/V  -  0,  the  orbit  is  a  nure  circle 
of  radius  r  passing  over  the  target.  Hence,  the  position  is 
evenly  distributed  with  time  about  the  circle,  and  the  value 
of  Rc  is  r  VT  or  R,./r  =  1.41.  Eq  (4°/  gives  the  value  Rr./r  =  1.35, 
a  difference  of  5%.  For  W/V  =■  1.0  tfhich  has  no  real  significance 
except  in  the  concept  of  a  mathemat i c al  limit,  Eq  (49)  gives 
Rp/r  =  3.18.  A  reasonable  assumption  for  the  limit  of  Rc  as 
W/V  approaches  1.0  is  yjj/2  .  The  value  of  y^  given  by  Eq  (16) 
of  Sectio  2a  is  -20tr  which  is  also  the  value  of  yj? .  Hence, 
theoretical  1  v  Rc  =*TTr  and  Rc/r  “TT=  3.14.  The  result  of  Eq  (49) 
differs  by  1%. 


T!:ii  s  ,  it  may  be  presumed  that  Eq  (49)  is  an  accurate 
method  for  determining  the  radius  of  the  CEP.  Then,  if  Rc> 
W  and  V  are  considered  design  requirements  for  the  system, 
the  solution  of  Eq  (48)  for  r, 

r  -  Rc  (1.27  -  W/V)  (50) 

1.79  -  0.84  W/V 

gives  the  necessary  turning  radius  for  the  system. 
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4 .  Concluslon3 


The  coordinates  of  four  characteristic  points  and  the 
total  downwind  displacement  defining  the  geometry  of  a  typical 
descending  orbit  can  be  calculated  directly  as  functions  of 
turning  radius  and  the  ratio  of  wind  velocity  to  horizontal 
system  velocity  relative  to  air  (Eq's  2, 3, 5, 7  and  14-18).  The 
total  time  elapsed  per  orbit,  is  equal  to  the  time  required 
to  turn  a  complete  circle  in  calm  eir  at  a  speed  equal  to  the 
minimum  ground  speed  in  a  wind  (Eq  Id) .  Points  at  equal  time 
intervals  along  the  ground  track  may  be  computed  (Eq's  35-38) 
with  the  time  interval  chosen  as  a  function  of  the  oscillation 
of  the  system.  Thus,  effects  of  the  real  system  may  be  approx¬ 
imated.  Analysis  of  the  deployment  envelope  shows  that,  in 
general,  the  median  altitude  (Eq  45)  of  the  system  at  the  first 
pass  over  the  target  is  sufficient  to  ensure  more  than  one  com¬ 
plete  orbit  prior  to  impact  for  missions  utilizing  the  offset 
delivery  capability.  Hence,  the  assumption  that  the  target  zone 
is  dependent  on  orbit  geometry  is  valid.  Impact  position  prob¬ 
ability  may,  therefore,  be  estimated  according  to  the  hypothesis 
that  impact  will  occur  at  equal  probabilitv  with  time  during  a 
typical  orbit. 

Computer  iterations  show  good  agreement  with  basic  geo¬ 
metric  relations,  but  also  show  that  the  magnitude  of  the 
time  increment  may  have  a  significant  effect  on  calculated 
behavior  near  the  target.  Downwind  drift  during  long  crosswind 
approaches  increases  linearly  with  offset  distance  and  is 
inversely  proportional  to  glide  ratio  while  somewhat  directly 
proportional  to  turning  radius. 

The  linearity  of  probability  values  determined  from 
computer  results  allows  accurate  empirical  formulation  of 
the  relation  of  probability  to  radius  from  target  ,  turning 
radius,  wind  velocity,  and  horizontal  system  velocity  (Eq  47). 
Hence,  reliable  equations  showing  the  functional  relationships 
between  radius  of  the  CEP,  turning  radius,  wind  velocity,  and 
horizontal  system  velocity  may  be  formulated  (Eq's  49  and  50) 
which  are  suitable  to  mission  and  design  analysis.  From  the 
empirical  relationship,  we  see  that  increased  glide  ratio,  or 
horizontal  system  velocity, has  the  effect  of  decreasing  the 
radius  of  the  CEP  while  also  decreasing  the  wind  disf-ortion 
during  the  approach. 


APPENDIX 


Glossary  of  Symbols  and  Abbreviations 


Positive  or  negative  unity  based  on  homing 
tes  t . 


b 

CEP 

CW 

DW 

h 

ho 
L  h 
L 


X  O 

i/d 


m 


P 


Value  of  P  at  R/r  =  0  as  determined  by  Method 
of  Least  Squares  . 

Abbreviation  for  circle  of  equal  probability 
with  center  at  the  target. 

Initial  position  crosswind  of  target. 

Initial  position  downwind  of  target. 

Median  altitude  at  first  pass  over  target 

Nominal  deployment  altitude 

Error  in  deployment  altitude. 

Net  downwind  displacement  during  orbit  measured 
from  the  ordinates  of  the  most  upwind  position 
and  the  most  downwind  position. 

Nominal  range  at  deployment 

Error  In  range  at  deployment. 

Glide  ratio  equal  tr>  ratio  of  range  to  altitude 
loss  and  also  equal  to  ratio  of  lift  to  drag 
and  horizontal  to  vertical  velocity  comnonents  . 

Slope  of  curve  for  P  vs  R/r  as  determined  by 
Method  of  Least  Squares. 

Probability  of  lmosct  within  circle  of  radius  R. 


R 

Rc 


Radius  of  circle  with  center  at  target,  also 
referred  to  as  distance  from  target. 

Radius  of  the  CEP. 


r  Turning  radius  of  the  system,  considered  a 

constant  by  design. 

S  Square  of  distance  from  target  to  end  of  arc  d. 

T  Time  elapsed  during  a  typical  orbit  In  a  wind. 
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To 


Time  elapsed  during  a  complete  360-degree  turn 
in  Calm  air. 


t  Elapsed  time  at  any  position  during  an  orbit. 

UW  Initial  position  upwind  of  target. 

V  Horizontal  component  of  sy  ttem  velocity  relative 

V<j  Vertical  component  of  system  velocity  relative 

to  air,  also  referred  to  as  rate  of  descent. 

Vx  Component  of  V  oarallel  to  x-axis. 

Vy  Component  of  V  parallel  to  y-axis  . 

W  Wind  velocity  directed  toward  the  target  along 

the  positive  y-axis. 

W  Random  gust  velocity. 

x  Displacement  along  the  x-axis  which  is  in  the 

crosswind  direction. 


y  Displacement  along  the  y-axis  which  is  in  the 

upwind  direction. 

0^  Angle  between  V  and  pure  croiswind  resultant 

of  V  +  W  at  most  upwind  position  on  orbit. 

9  Angle  of  arc  during  a  turn. 

0  Inclination  of  V  from  the  x-axis. 

Increment  of  tin-'  o'msed  during  turn  through 
arc  9,  also  referre'1  to  as  time  lag. 

(  )'  Superscript  denoting  value  at  end  of  arc  9. 

(  or  (  >2  Subscripts  denoting  value  at  one  of  two  possible 

end  points  of  arc  9  resulting  from  homing  control. 

(  )w  Average  at  various  initial  positions. 

(  )ry-  Average  for  various  values  of  time  increment. 
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Fig  3.  Flow  Chart  for  Computer  Program 
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Fig  11  EXPANDED  PORTION  OF  CASE  6  SHOWING  INSTANTANEOUS  SYSTEM  ORIENTATIONS 
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Fig  27  CHARACTERISTIC  POINTS  ON  THEORETICAL  ORBITS  CORRESPONDING  TO  CASES 
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Fig  29  LEAST  SQUARES  ANALYSIS  OF  P  vs  R/r  AT  VARIOUS  W/V 
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